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Abstract 



We discuss magnetic monopole solutions of the Einstein- Yang-Mills-Higgs equations with a positive 
cosmological constant. These configurations approach asymptotically the de Sitter spacetime background 
| and exist only for a nonzero Higgs potential. We find that the total mass of the solutions within the 

cosmological horizon is finite. However, their mass evaluated by using the surface counterterm method 
QO ' outside the cosmological horizon at early/late time infinity generically diverges. Magnetic monopole 

£ — . solutions with finite mass and noninteger charge exist however in a truncation of the theory with a 

vanishing Higgs field. Both solutions with a regular origin and cosmological black holes are studied, 
C****" ' special attention being paid to the computation of the global charges. 

O ' 

1 Introduction 

o 

As found some time ago by 't Hooft and Polyakov, certain gauge theories admit classical, particle-like 
b , solutions with quantised charge and finite energy pQ. The magnetic monopoles inevitably arise in grand 

unification theories whenever the spontaneous symmetry breakdown of this theory generates a U(l) subgroup 
and are stabilised by a quantum number of topological origin, corresponding to their magnetic charge. Since 
the Standard Model of particle physics contains a U(l) subgroup, monopoles with masses of the order 10 17 
GeV are predicted which leads to the so-called "monopole problem" . The accepted solution to this problem 
■ is the theory of inflation. 

When coupling the Yang-Mills-Higgs (YMH) system to gravity, a branch of globally regular monopoles 
emerges smoothly from the corresponding flat space solutions. The nonabelian black hole solutions emerge 
from the globally regular configurations, when a finite regular event horizon radius is imposed |21E]- These 
solutions cease to exist beyond some maximal value of the coupling constant a, which is proportional to the 
ratio of the vector meson mass and Planck mass. 

It has been speculated that such configurations might have played an important role in the early stages of 
the evolution of the Universe. Also, various analyses indicate that the monopole solutions are important in 
quantum theories. However, most of the investigations in the literature correspond to an asymptotically flat 
(AF) spacetime. Because of the physical importance of these objects, it is worthwhile to study generalisations 
in a different cosmological background. 

As discussed in [3], the properties of magnetic monopole solutions in Anti-de Sitter (AdS) spacetime 
are rather similar to their asymptotically flat counterparts. Nontrivial solutions exist for any value of the 
cosmological constant A < 0; as a new feature, one finds a complicated power decay of the fields at infinity 
and a decrease of the maximally allowed vacuum expectation value of the Higgs field. 

For a positive cosmological constant, the natural ground state of the theory corresponds to de Sitter (dS) 
spacetime. Solutions of Einstein equations with this type of asymptotics enjoyed recently a huge interest in 
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theoretical physics for a variety of reasons. First at all, the observational evidence accumulated in the last 
years (see, e.g., ref. |S|), seems to favour the idea that the physical universe has an accelerated expansion. 
The most common explanation is that the expansion is driven by a small positive vacuum energy (i.e. a 
cosmological constant A > 0), implying the spacetime is asymptotically dS. Furthermore, dS spacetime plays 
a central role in the theory of inflation (the very rapid accelerated expansion in the early universe), which is 
supposed to solve the cosmological flatness and horizon problems. 

Another motivation for studying solutions with this type of asymptotics is connected with the pro- 
posed holographic duality between quantum gravity in dS spacetime and a conformal field theory (CFT) 
on the boundary of dS spacetime [Hj. The results in the literature suggest that the conjectured dS/CFT 
correspondence has a number of similarities with the AdS/CFT correspondence, although many details 
and interpretations remain to be clarified (see |S] for a recent review and a large set of references on this 
problems). 

In view of these developments, an examination of the classical solutions of gravitating fields in asymp- 
totically dS spacetimes seems appropriate, the physical relevant case of a spontaneously broken nonabelian 
gauge theory being particularly interesting. 

Considering a static coordinate system (which generalises for A > the usual Schwarzschild coordinates) , 
these solutions will present a cosmological horizon for a finite value of the radial coordinate, where all 
curvature invariants stay finite. Similar to the well-known (electro)-vacuum solutions, there is no global 
timelike Killing vector, as the norm of the Killing vector d/dt changes sign as it crosses the cosmological 
horizon. Inside the horizon, the Killing vector is timelike, and this can be used to calculate the conserved 
charges and action/entropy inside the cosmological horizon Outside the cosmological horizon, however, 
with the change of d/dt to a spacelike Killing vector, the physical meaning of the conserved quantities is 
less clear. 

Recently a method for computing conserved charges (and associated boundary stress tensors) of asymp- 
totically dS spacetimes from data at early or late time infinity was proposed |l(Jj . in analogy with the 
prescription used in asymptotically AdS spacetimes This approach uses counterterms on spatial bound- 
aries at early and late times, yielding a finite action for asymptotically dS spacetimes. The boundary stress 
tensor on the spacetime boundary can also be calculated, and a conserved charge (spacelike, due to its 
association with the Killing vector d/dt) - now interpreted as the mass of the solutions, can be defined. 

The main goal of this paper is to present a study of the basic properties of the spherically symmetric 
monopole solutions in Einstein- Yang-Mills-Higgs (EYMH) theory with a positive cosmological constant. 
(Solutions of the Abelian-Higgs model with A > have been studied in 52 •) Although an analytic solution 
to the EYMH equations appears to be intractable, we present both analytical and numerical arguments for 
the existence of both solutions with a regular origin and cosmological black-hole, approaching asymptotically 
the dS background. As found in ^3], the features of the dS solutions of a spontaneously broken nonabelian 
gauge theory are rather different as compared to the asymptotically flat or AdS counterparts. The most 
interesting feature is that the mass of dS magnetic monopoles evaluated at timelike infinity according to 
the counterterm prescription generically diverges, although the total mass within the cosmological horizon 
of these configurations is finite. Also, no solutions exist in the absence of a Higgs potential. 

However, in the last years it became clear that a nonasymptotically flat metric background may allow 
for nonabelian magnetic monopole solutions even in the absence of a Higgs field. For example, a charge- 
one monopole solution has been found by Chamseddine and Volkov 14 in the context of the Af = 4 
D = 4 Frecdman-Schwarz gauged supergravity [15| . This is one of the few analytically known configurations 
involving both non-abelian gauge fields and gravity, the expression of the magnetic potential coinciding, in a 
suitable coordinate system, to that of the Prasad-Sommerfield monopole [17| . A solution with many similar 
properties exists also JS| m a version of Af = 4 d = 5 Romans' gauged supergravity model with a 
Liouville dilaton potential. 

As found in |2()l 121) . monopole- type solutions exist even in a simple Einstein- Yang-Mills (EYM) theory 
with a negative cosmological constant. These asymptotically AdS nonabelian solutions have finite mass and 
an arbitrary value of the gauge potential at infinity, i.e. an arbitrary value of the magnetic charge. There 
are also finite energy solutions for several intervals of the shooting parameter (the value of gauge function 
at the origin or at the event horizon), rather than discrete values and stable monopole solutions in which 
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the gauge field has no zeros. 

It is natural to conjecture the existence of similar configurations for A > 0. Solutions with a regular 
origin of the EYM-SU(2) system with positive cosmological constant have been considered by several authors 
(see [H]-HS], and also the systematic approach in Similar to the AdS case, the asymptotic value of 

the gauge potential for solutions with A > is not fixed, implying the existence of a nonvanishing magnetic 
charge. However, all asymptotically dS configurations are unstable. In this paper we propose to reconsider 
the properties of these EYM-A solutions, viewing them as magnetic monopoles in a model without a Higgs 
field. Both solutions with a regular origin and cosmological black hole configurations are discussed. 

The mass and boundary stress-tensor as well as the thermodynamic quantities of both EYMH and EYM 
solutions discussed in this paper are computed by using the counterterm formalism proposed in |10j . We 
also comment on the implications of these solutions for the conjectured dS/CFT correspondence. 

The paper is structured as follows: in the next Section we present the general framework and analyse 
the field equations and boundary conditions. In Section 3 we present our numerical results. In Section 4 we 
compute the solutions' mass and action, while in Section 5 we comment on the general picture in a different 
coordinate system. We conclude with Section 6, where our results are summarised. The Appendix generalises 
the construction of [2Z| for a positive effective cosmological constant. We show that the BPS monopole 
of four-dimensional YMH theory continue to be exact solutions even after the inclusion of gravitational, 
electromagnetic and dilatonic interactions, provided that certain non-minimal interactions are included. 



2 General framework 
2.1 Action principle 

The action for a gravitating non-Abclian SU (2) gauge field coupled to a triplet Higgs field is 

iB = L d4xV ~ 9 - 2A ) + c -) - L d " x ^ w 



with Newton's constant G and 



C M = -^Tv(F^F^) - ^Tt(D^D^) - (2) 



where V($) = ±ATr(<I> 2 — T] 2 ) 2 is the usual Higgs potential, with r\ the vacuum expectation value (vev). The 
gauge field strength tensor is given by 

F liU = d^A u -d„A li + ie[A ll ,A u \, (3) 

with Dfj, = d^ + ielAft, ] being the covariant derivative, A^ the gauge potential and e the Yang-Mills coupling 
constant. 

The last term in 0) is the Hawking-Gibbons surface term [2Hj, where K is the trace of the extrinsic 
curvature for the boundary dM. and h is the induced metric of the boundary. Of course, this term does not 
affect the equations of motion but is important when computing the mass and action of solutions. 

Apart from the cosmological constant, the theory contains three mass scales: the Planck mass Mpi = 
1 / \/G, the gauge boson mass M\y = ev and the mass Mr = y/Xr] of the Higgs field. 

Varying the action (JJJ with respect to g^ v , A^ and <I> we have the field equations 

-Li) M (v^H = \ie[W*}> (4) 



^D„{y/=gD»*) + A($ 2 - r) 2 )$ = 0, 
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where the stress-energy tensor is 

TjS ] = Tr^D^D^-^D^D^-g^V^), (5) 
= 2Tr(F m F u0 g^ - \g, v F a pF^) . (6) 
The EYM model is found in the limit of a vanishing Higgs field (f> = 0, V(Q) = 0. 
2.2 The ansatz 

For A < 0, the above field equations present a large variety of solutions, including configurations with axial 
symmetry only |29j. |30| . However, here we will restrict to the spherically symmetric case. 

When discussing physics in dS spacetime, there are many coordinate systems to choose among. In this 
work will consider mainly the static system. The advantage of these coordinates is their obvious simplicity 
and the time independence. This coordinate system is also computationally more convenient, since we will 
deal with ordinary differential equations. The disadvantage is that the expansion of the universe is not 
manifest. (In Section 5 we will discuss our results in a cosmological coordinate system. Both can be used to 
approach the timelike infinity and will lead to different form of the boundary metric.) 

Therefore we consider a line element given by 

df 

ds 2 = —-- + r 2 (d0 2 + sin 2 edifi 2 ) - a 2 (r)N(r)dt 2 (7) 
N(r) 

where a convenient parametrisation of the metric function N(r) is 

„ T , , „ 2m(r) Ar 2 

N(r) = 1 ^ - — . 8 

r S 

The function m(r) is usually interpreted as the total mass-energy within the radius r. For m(r) = 0, the 
empty dS space written in a static coordinate system with a cosmological horizon at r = r c = y/S/A. is 
recovered. 

For the SU(2) Yang-Mills field we use the minimal spherically symmetric (purely magnetic) ansatz em- 
ployed also in previous studies on asymptotically flat (or AdS) solutions, in terms of a single magnetic 
potential w(r) 



A 

2e 



Y |w(r)ndfl + ( cos 6»r 3 + u(r)T 2 sin 9) dtp}, (9) 
while the Higgs field is given by the usual form 

$ = <t> T 3) (10) 

where n are the Pauli matrices. 



2.3 The equations of motion 

For a nonvanishing asymptotic magnitude of the Higgs field 77, dimcnsionlcss quantities are obtained by 
considering the rescaling 1 

r—*r/(rje), m— >m/(r]e), <3> — > 0?7, A— >A(e??) 2 . (11) 

As a result, the EYMH equations then depend only on the dimcnsionlcss parameters a = Mw /(eMpi), 
(5 = Mh /Mw and on A. 

It can be proven that a vanishing asymptotic magnitude of the Higgs field 77 implies $ = 0. 
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Within the above ansatz, the reduced EYMH action can be expressed as 



dr dt 



am 1 - a 2 a(u' 2 N + ^ ^ + -<j> ,2 r 2 N + u 2 <j> 2 + VU)r 2 ) , (12) 
2r 2 2 



where the prime denotes the derivative with respect to the radial coordinate r and V(<p) = (3 2 ((f> 2 — l) 2 /4, 
with P 2 = \/e 2 . 

The EYMH equations reduce to the following system of four non-linear differential equations 



(AW)' = acu( {uj2 1} + cb 2 ), (Nar 2 ^')' = a{2u?<\> + r 2 ^-). (14) 



In the absence of the Higgs field, one defines a set of dimensionless variables by performing the following 
rescalings r — > (V4~7rG/e)r, A — * (e 2 /47rG)A and m — > (eG/\/47rG)m. The equations of the EYM-A system 
model then read 

m' = cu' 2 N+ {Uj2 ~ 1) \ a' = --au,' 2 , ( W)' = ~ 1} . (15) 

2r r r 

2.4 Known solutions 

Several explicit solutions of the above equations are well known. The Schwarzschild-de Sitter solution 
corresponds to 

2M Ar 2 

w(r)=±l, o-(r) = l, 0(r) = with u = 0, AT(r) = 1 , (16) 

r 3 

and describes a black hole inside a cosmological horizon as long as N(r) has two positive zeros, i.e. M < 
l/3\/A. The Nariai solution is found for M = 1/3 \/A, in which case the function N(r) has a double zero at 
r c = l/VX, corresponding to the position of an extremal cosmological horizon. 
The second solution corresponds to an embedded U(l) configuration with 



w(r) = 0, o-(r) = 1, 0(r) = 1, ^V(r) = 1 + — - — (17) 



2M . a 2 Ar 2 

T 



and describes a particular parametrisation of the (magnetic-) Reissner-Nordstrom-dS (RNdS) black hole. 
Since this solution proves to be important in understanding the features of the nonabelian configurations, 
we present here a brief review of its properties (see e.g. [HJj |32| . pTH] for more details). 

For A > 0, up to four zeros of N(r) can exist. Three of the four zeros correspond to horizons since 
the first zero has always negative value and thus has no physical meaning. The two inner horizons r_, r + 
with r_ < r + correspond to the well known Cauchy, respectively event horizon of the Reissner-Nordstrom 
solution, while the third outer horizon r c > r + is the cosmological horizon. Extremal black hole solutions - 
like in the AF case - are also possible. Then, we have r_ = r + = with N{rh) = jV'| r=rh = 0. This leads 
to the equation: 

Ar 4 h -r 2 +a 2 =0, (18) 

which is solved by 



ru/c= vk^ 1±Vl_4a2A for i- A>0 - (19) 



Obviously, the appearance of horizons in dS space is restricted by a 2 < i. The corresponding mass of the 
extremal solution is given by: 

M e 2^ + I£. (20 ) 
ir h 3 
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As the parameter M increases (relative to a, with M > 0), the outer black hole and cosmological horizons 
move closer together. The charged Nariai solution is obtained when these horizons coincide at = r c ; this 
is the largest charged asymptotically dS black hole. For a given a, the charged Nariai black hole has the 
maximal mass. 

For completeness, we present also another solution which is less relevant to the situation discussed in 
this paper. This solution is found for a particular value of the cosmological constant A = 3/(167rG) and a 
vanishing Higgs field <f> = (3 = 0. It reads 

ds 2 = i * r * Ar2 +rW-dt 2 , u{r) = ^l- 2 -Kr\ (21) 

corresponding to an S 3 x R Einstein universe. 



2.5 Boundary conditions for the EYMH model 

We want the generic line element (JJJ to describe a nonsingular, asymptotically de Sitter spacetime outside 
a cosmological horizon located at r — r c > 0. Here N(r c ) — is only a coordinate singularity where all 
curvature invariants are finite. Outside the cosmological horizon r and t changes the character (i.e. r becomes 
a timelike coordinate for r > r c ). A nonsingular extension across this null surface can be found just as at 
the event horizon of a black hole. The regularity assumption implies that all curvature invariants at r = r c 
are finite. Also, all matter functions and their first derivatives extend smoothly through the cosmological 
horizon, e.g. in a similar way as the U(l) electric potential At — Aq + Q/r of a RNdS solution. 

Similar to the case A < 0, it is natural to consider two types of configurations, corresponding in the usual 
terminology to particle-like and black hole solutions. 

For particle-like solutions, r = is a regular origin, and we find the following behaviour: 

m(r) = ~a 2 (12b 2 + 3a 2 + 2/3 2 )r 3 + 0(r 5 ), a(r) = a + \(a 2 + 8b 2 )a a 2 r 2 + (9(r 3 ), 
6 2 

to{r) = 1-br 2 + 0(r 4 ), <j>{r) = ar + (9(r 3 ), (22) 
where a and b are free parameters which are found by solving the field equations. 



2.5.1 Expansion at the horizon 

The cosmological horizon is located at some finite value of r, where N(r c ) = 0, <r(r c ) ^ 0, the gauge potential 
and the scalar field taking some constant values. The corresponding expansion as r — > r c is 

m(r) = y(l-^) + mi(r-r c ) + 0(r-r c ) 2 , a(r) = a c + a x (r - r c ) + 0{r - r c ) 2 , 

u(r) = uj c + u) X {r-r c ) + 0(r-r c f, <j)(r) = C + <fn(r - r c ) + 0(r - r c f , (23) 

where 



;-^K 2 - 1) 2 + ^ 2 C + v{d> c y c ), co, = ± C c{u) \ l) 



1 /„ 9 ,9 idV 



Niri 



^ a c 2 1 2 | 1 ±2 2\ 



(2^+r 2 — ), a x = -at'{wi + -<fiii), (24) 

r c 2 



a c , lu c , 4> c being arbitrary parameters and iVi = N'(r c ) = (1 — 2mi — Ar 2 )/r c . Note that the condition 
N'(r c ) < should be satisfied, which imposes the following constraint 

^ 2 C - I) 2 + ulcfc + V^ c )r 2 c > i-M. (25) 



G 



We will consider also cosmological black hole solutions. These configurations possess an event horizon 
located at some intermediate value of the radial coordinate < < r c , all curvature invariants being finite 
as r — > rh- As r — > rh, the functions m, a, ui and <p present an expansion very similar to l|23|) 



m(r) 
uj(r) 



(1 



Ar 



) + rn'(r h ) + 0{r - r h ) 2 , a{r) = a h + a'(r h )(r - r h ) + 0(r - r h f 



= u; h + Lu'(r h ){r-r h ) + 0(r-r h ) 2 , 0(r) = ^> h + <jt \r h ){r - r h ) + 0(r - r h f 



(26) 



where 



1 



1 



r L) h (uJ 2 h - 1) 



m ' { r h) = « 2 (^|K - I) 2 + + V(<t> h )ri), J{r h ) = ^--i ~' + Uh <ft), 

N'ir h )r s h d<p r h r h 2 



(27) 



ah, 4>h being arbitrary parameters and N'(rh) = (1 — 2m! {rh) — Ar^/r^. The obvious condition 



N'(rh) > imposes the constraint 

1 



l-Ar 2 



2r2 {ul-l? + u 2 h <i>l+V{<t> h )rl< 2 



(28) 



Both the event and the cosmological horizon have their own surface gravity k given by 

<c = -\9 U 9 rr {dr9uf 



the associated Hawking temperature being Th — l K l/(27r)- 



2.5.2 Expansion at infinity 

The analysis of the field equations as r — > oo implies a more complicated picture as compared to the AF or 
asymptotically AdS case, since the cosmological constant enters in a nontrivial way the solutions' expression 
at infinity. 

We suppose that the Higgs scalar approaches asymptotically its vev, while the magnetic gauge potential 
w(r) vanishes. This assures the absence of supplementary contributions to the cosmological vacuum energy 
(apart from A) as r — > oo, while the gauge field approaches asymptotically the U(l)-Dirac monopole field. 
The asymptotic expression of the metric functions m(r), a(r) and the matter functions w(r) and (j)(r) is 
found by finding an approximate solution of the field equations <|13ll - H14[) for these boundary conditions. 

Here it is instructive to work with dimensionfull variables, the corresponding expression after considering 
the rescaling being straightforward. 

A systematic analysis of the matter field equations reveals that a positive cosmological constant sets a 
mass bound for the gauge sector, Mf, = y/A/12. The expression of the gauge field as r — > oo for Mw < Mb, 
which leads to finite mass solutions is 

w(r)~ Cl r k \ with k l = ~(l + ^^M 2 ^\ . (29) 

This contrasts with the exponential decay found in an asymptotically flat spacetimc. 
For Mw > Mb, the large r behaviour of the gauge field is 

uj(r) - c 3 r~^ 2 sin^x (r) , with 9i{r) = (Pi logr + c 4 ), Pi - ^M w /M 2 - 1, (30) 

which we will find leads to a logarithmic divergence in the asymptotic expression of the mass function m(r). 

The analysis of the scalar field asymptotics is standard. Strominger's mass bound jH| (i-e. the dS version 
of the Breitenlohner-Freedman bound 0) is M| = 3 A/4 and separates the infinite energy solutions from 
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solutions which may present a finite mass (this would depend also on the gauge field behaviour). For small 
enough values of the Higgs field mass, Mh < Ms the scalar field decay leading to a finite asymptotic value 
of m(r) is 



(r) ~ rj + c 2 r k2 , with k 2 = ~ (l + y/l - M 2 H /M 2 



(31) 



We found numerical evidence for the existence of a secondary branch of solutions decaying as 



(r) ~ ri + c 2 r' k2 , with % = -\ (l- ^l- M%/M^ . (32) 



However, this decay leads to a divergent value of the function m(r) as r — > oo. 

For a Higgs mass exceeding Strominger's bound, the scalar field behaves asymptotically as 



cj)(r) ~ n + c 5 r- i/z sin tt 2 (r) , with # a (r) = (P 2 logr + ce), P 2 = -^M^/Afj - 1, (33) 

(the constants Cj which enter the above relations are free parameters). 

One may think that this bound may be circumvented by solutions with a vanishing Higgs potential. 
However, by rewriting the Higgs field equation in the form 

-(Nar 2 U 2 )')' = o(Nr 2 ^ 2 + 2w 2 tf + r 2 <A, 
2 dcp 

and integrating it between the origin and the cosmological horizon, it can easily be proven that no nontrivial 
solutions exist for V(4>) = or for a convex potential. 

Once we know the asymptotics of the matter fields, the corresponding expression for the metric functions 
results straightforwardly from the equations l|13l) . 

For Mh < Mg and Myj < Mf, the mass of solutions tends to a constant value as r — > oo and we find the 
asymptotic expressions 

- . , 2 A c 2 r 2kl+1 M 2 Ak 2 c 2 r 2k *+ 3 1 . 



fc 2 c 2 r 2fc i~ 2 , k 2 c 2 2k2 
2g 2 (k 1 -l) 



loga(r) 8ttG(^ _ + 212^) + . . . . (35 ) 



In the case My/ > Mb the metric function m(r) gets a divergent contribution from the nonabelian field, 
whose leading order expansion is 

ttGc 2 ( A 
m YM (r) ~ AttGMx -2 (-2*1 + sin2tfi)n 2 + — ^(8Pi cos 2 *i (36) 

Pi V I2cr 

+2(1 + 4P! 2 )*! + (4P X 2 - 1) sin 2*i 

(note the presence in this expression of a divergent term originating in the kinetic Higgs term). 

For a Higgs field mass exceeding the Strominger bound Mh > Ms, we find a very similar asymptotic 
form of the mass function m(r), presenting the same type of divergencies 

m H (r) ~ 4nGM 2 - < ^f- (\t?{-2^ 2 + sin2* 2 ) + ^(24P 2 cos 2 * 2 (37) 
4P 2 \ 6 



+2(9 + 4P|)tf 2 + (-9 + 4P 2 ^) sin 2* 

Since the order of the cosmological constant is believed to be much smaller than both the gauge and Higgs 
boson mass, our results indicate that monopoles in the universe will always have a divergent mass function 



8 



m(r) as evaluated at timelike infinity. Moreover, as discussed in Section 4, we could not find numerical 
solutions with Mh < Ms and iV% < Mf,. 

One should also study the asymptotic expression of the a function. Supposing Mh > Ms and My/ > Mb 
one finds 



loga(r) ~ ^ 



c 



.3 

3 ((l + 4 J P 1 2 )(9 + 4P 1 2 ) + 3(-3 + 4P 1 2 )cos2* 1 -6Pi(5 + 4P 1 2 )sin2*i) (38) 



9 + 4P 2 

+c§(9 + 4P| - 9 cos 2* 2 - 6P 2 sin 2* 2 



The corresponding expressions for Mh > Ms or Mw > Mb can easily be read from lf3%)l. (|3*5|l . 

The solutions with My = Ms, = Mft saturate these bounds and lead also to infinite mass configu- 
rations. The EYMH equations lead to a matter fields expression 

W (r)~^ + 5^I, + (39) 



(with C7, .., cio real constants), while for metric functions one can write m(r) = m (r) + m (r), where 

/A c 3 \ 

m YM (r) = 4^Gf-^-^(c 7 -2c 8 + c 8 logr) 3 +r; 2 (c?logr + c 7 c 8 logr + ^log 3 r)J +... (40) 

m ff (r) = AttG (-—^—(-2c w + 3c 9 + 3c w logr) 3 + ^yy 2 log r (3c 2 + 3c 9 ci logr + c w log 2 r) 

\ ilDCio D 



and 



2irG 



logo-(r) = -7^3 (9c 7 -30c 7 c 8 + 26c 2 + 3c 8 log r(6c 7 - 10c 8 + 3c 8 logr)) (41) 

7rG 

- 7^-3 (2ci - 6c 9 cio + 9cg + 3ci log r (-2c w + 6c 9 + 3ci log r)) + . . . . 

One can see that even for a logarithmic diverging mass function, the spacetime is still asymptotically dS, 
presenting the same conformal structure at infinity as the vacuum dS solution, since gtt ~ —1 + Ar 2 /3 + 
0(\ogr /r). Also, the Ricci scalar stays finite as r — > 00, R — > 4A. 

Solutions with A < approaching at infinity the AdS background despite the presence of a divergent 
ADM mass have been considered by various authors in the last years, mostly for a scalar field matter content. 
Restricting to the case of models with nonabelian matter fields in the bulk, we mention the SU(2) hairy 
black holes in [33], the family of globally regular solutions in D=4, JV = 4 gauged supergravity in |33], and 
the self-gravitating Yang-monopoles in [3j)|. The situation for A > is less studied. However, it is natural 
to expect the existence of solutions with similar features for a dS background too, the EYMH monopoles 
being a rather complicate case. 

2.6 Boundary conditions for the EYM model 

The field equations imply the following behaviour for r — > in terms of two real parameters b, a 

w(r) = l-6r 2 + 0(r 4 ), m(r) = 2b 2 r 3 + 0(r 4 ), a(r) = cr + <ib 2 (T r 2 + <3(r 4 ). (42) 
The corresponding expansion near the comological horizon is 

m(r ) = f (l- ^) + ^f(r ~ r„), «(r) = u„ + (1 _™$'$_ 1)a (r - r c ), (43) 
<j(r) =a c - *2L u \r h ) 2 (r - r c ) + 0(r - r h ) 2 , 
where w c , a c are real parameters. 
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When discussing the pure EYM system with A > 0, there are no restrictions on the asymptotic value of 
the gauge potential |22) . The field equations imply the following expansion at large r 

m(r) =M+(*£- \H - If) \ + O (i) , w(r) = .0 + £ + O (1) , 

ff (r) = l-^+o(l), (44) 

where -^f and Ci are constants determined by numerical calculations. 

For cosmological black hole solutions having a regular event horizon at r = rh > 0, we find the following 
expansion near the event horizon 

m(r) = ^(l-^)+^4^(r-r fc ), u{r) = u h + (1 ^g r l"jg ) _ 1) , (r - r h ), (45) 
= <Th - T^oj'(r h ) 2 (r - r h ) + 0(r - r h f, 

with Wh, Oh real parameters. 



3 Counterterm method and conserved charges 
3.1 General formalism 

The computation of the conserved charges including mass in an asymptotically dS spacetime is a difficult 
task. This is due to the absence of the spatial infinity and of a globally timclikc Killing vector in such a 
spacetime. One prescription to compute conserved charges in asymptotically dS was developed by Abbott 
and Deser |37l . In this perturbative approach one considers the deviation of metric from pure dS space which 
is the vacuum of the theory and measure the energy of fluctuations to find the mass. 

In |10| . a novel prescription was proposed, the obstacles mentioned above being avoided by computing 
the quasilocal tensor of Brown and York (augmented by the AdS/CFT inspired counterterms) , on the 
Euclidean surfaces at future/past timelike infinity 1 ± . The conserved charge associated with the Killing- 
vector d/dt - now spacelike outside the cosmological horizon- was interpreted as the conserved mass. This 
allows also a discussion of the thermodynamics of the asymptotically dS solutions outside the event horizon, 
the boundary counterterms regularising the (tree-level) gravitation action. The efficiency of this approach 
has been demonstrated in a broad range of examples. 

A thorough discussion of the general formalism has been given e.g. in |38| . and so we only recapitulate 
it here. In this approach one starts by considering the general path integral 

(<72,* 2 ,S2|fli,*i,Si) = [ D\g,y}exp(iI[g,V}), (46) 

which represents the amplitude to go from a state with metric and matter fields [gi, on a surface Si to 
a state with metric and matter fields [g2, ^2] on a surface S2. The quantity D [g, 'J] is a measure on the 
space of all field configurations and I [g, ^] is the action taken over all fields having the given values on the 
surfaces Si and 82- For asymptotically dS spacetimes we replace the surfaces Si, £2 with histories Hi,H 2 
that have spacelike unit normals and are surfaces that form the timelike boundaries of a given spatial region. 
The amplitude l|4^|) describes quantum correlations between differing histories [gi, and [<?2, ^2] of metrics 
and matter fields, with the modulus squared of the amplitude yielding the correlation between two histories. 

In this approach, the initial action JQ) is supplemented by the boundary counterterm action I ct depending 
only on geometric invariants of these spacelike surfaces. I c t regularizes the gravitational action and the 
boundary stress tensor. In four dimensions, the counterterm expression is (in this section we do not consider 
rescaled quantities; also to conform with standard conventions in the literature on this subject, we note 
A = 3/£ 2 ) 
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with R the curvature of the induced metric hij and Jg M ± indicates the sum of the integral over the early 
and late time boundaries. In what follows, to simplify the picture, we will consider the X + boundary only, 
dropping the ± indices (similar results hold for T~ ) . 

The boundary metric can be written, at least locally, in a ADM-like general form 

ds 2 = h l3 dx l dx 3 = Nfdt 2 + a ab {dijj a + N a dt) (dvj/ + N b dt) , (48) 

where Nt and N a are the lapse function and the shift vector respectively and the ip a are the intrinsic 
coordinates on the closed surfaces E. Varying the action with respect to the boundary metric hij, we find 
the boundary stress-energy tensor for gravity 

T «-^- (49) 



the corresponding expression for four dimensions being 39 



3 8ttG V 3 3 I 



( 2 /.. 1 
2 



(50) 



where Kij and are the extrinsic curvature and the Ricci tensor of the boundary metric, respectively. 
In this approach, the conserved quantity associated with a Killing vector on the X + boundary is given 

by 

05 = £d"^n%e, (51) 

where n % is an outward-pointing unit vector, normal to surfaces of constant r. Physically, this means that a 
collection of observers, on the hypersurface with the induced metric hij, would all measure the same value 
of 0£ provided this surface has an isometry generated by £\ 

If d/dt is a Killing vector on E, then the conserved mass is defined to be the conserved quantity 371 
associated with it. 

A tree-level evaluation of the path integral (|46|l for the EYM(-H) system may be carried out along the 
lines described in ref. for the vacuum case (see also the discussion in for the case of a gravitating 
U(l) field). Since the action is in general negative definite near past and future infinity (outside of a 
cosmological horizon) , we analytically continue the coordinate orthogonal to the histories to complex values 
by an anticlockwise 7r/2-rotation of the axis normal to them. This generally imposes, from the regularity 
conditions, a periodicity (3 of this coordinate, which is the analogue of the Hawking temperature outside the 
cosmological horizon. 

This renders the action pure imaginary, yielding a convergent path integral 

Z'= f e +1 (52) 

since I < 0. In the semi-classical approximation this will lead to h\Z' = +I c i. 
For a canonical ensemble with fixed temperature we can write 

F = Wl-TS, (53) 

where F is the canonical potential and T = 1/(3. For a converging partition function, we have F = I c i/j3 
and thus we find for the entropy of the system 

S = pTl-I cl . (54) 
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3.2 Spherically symmetric EYM(-H) solutions 
3.2.1 The boundary stress tensor 

The application of the general formalism to spherically symmetric EYM(-H) solutions discussed in the 
previous sections is straightforward. Working outside the cosmological horizon, we set following [3^1 r = r 
and rewrite the metric ansatz (JJJ as 

ds 2 = -/(r)rfr 2 + ^frdt 2 + rW , (55) 
J\ T ) 

where 

2m(r) 



f(T)=[p+—^-l) ■ (56) 



We choose dM. to be a three surface of fixed t > r c , which gives rii = 1/ ' yf{T)S T i). The extrinsic curvature 
Kij = h^Vkfij has the nonvanishing components 

t rsin 2 9 o(of - 2 fa') 

Kee = -7WY W = ~7WV u = v*<* ' (57) 

the corresponding expression for the trace of the extrinsic curvature being (in this section we do not consider 
rescaled quantities; also the prime here denotes the derivative with respect to r) 

2 of - 2 fa' , . 

K = + 9f3 /2 ■ 58 

to 2f 6 / z o 

For EYM solutions or EYMH configurations with M H < M s and M w < M b we find from JSUJ) that the 
nonvanishing components of the boundary stress-tensor are 



1 IM ( 1 \ , 1 Ml ( 1 



The mass of these solutions measured at the far future boundary of dS space, as computed from l|51|) , is 

Tt = —M, (59) 

where M is the parameter entering the asymptotic expansions 1)340. (|44f) . 

As discussed in the next section, all EYM solutions we found have M > (both black holes and particle 
like solutions). (Note that we didn't find any EYMH solutions with finite mass.) Thus 971 is negative, 
consistent with the expectation ^D] that pure dS spacetime has the largest mass for a singularity-free 
spacetime 071 = 0. If there is a CFT dual to a magnetic monopole, this mass translates into the energy of 
the dual CFT living on a Euclidean cylinder Rx S 2 . 

The boundary stress-tensor of the EYMH solutions with Mr > Ms or Mw > Mb (the only configurations 
we could find numerically is) 

T o^ T ^A TZ 9 1 iA + (^z\ T t = 1 9 1 Q_ +0 (logT 



" v 8ttG r 3 V t 4 J' 4 4ttG r 3 

where <7i(t) can be read from l|37() . I|36|l presenting a logarithmic divergence in r. This implies a divergent 
mass as computed from lf5l|) . 
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3.2.2 The magnetic charge 

The computation of the magnetic charge in the absence of a spacelike infinity is also difficult. However, 
as similar problems appear already in the abelian case, we may use the solution proposed in |33|. which 
generalises the methods of ref. 0U| to dS case. Working again outside the cosmological horizon, one defines 
the usual 't Hooft field strength tensor 

= Tr{$F^ - ±D„*D V *}, (60) 

where $ = 3>/|3>|. The induced metric h^ u — g^ u +u fJi u I/ projects the electromagnetic field on a specific slice 
of the foliation, where the vector u v = — ^= = 8" is normal to the induced metric h^u and n % = y/\f (t) \8\ 

is the (timelike) unit vector normal to hypersurface E. The magnetic field with respect to a slice r = const. 
is Bi = ^eijk^Fjk an d the charge density at future infinity is 

Pq = VvrfBi, 

Integrating the charge density over the hypersurface S, we obtain the total magnetic charge at future/past 
infinity 

Qm = ^<j> dS^e^ a Tr{$F va }. (61) 

For the particular case of spherically symmetric solutions one finds Q.,„ = 1 2 . The interpretation of the 
magnetic charge is analogous to that for 9JI noted above. Q m is the charge measured by a collection of 
observers following a given history; its conservation implies that observers following a different history would 
measure the same value of the magnetic charge 3 . 

The magnetic charge of the EYM solutions is evaluated at future/past infinity according to (see e.g.|21p 

1 f ,n ^-fnftt ^ T 3 



Qm = ^/ dS k V=gF kt = Q m ^, (62) 



where Q m = (1 — 



3.2.3 The entropy and action 

It can also be proven that the entropy of the monopole solutions associated with the cosmological event 
horizon is one quarter of the cosmological event horizon expected. 

Here we remark that, similar to the abelian theory |42j . the partition function is evaluated taking the 
magnetic charges as a boundary condition. Fixing the gauge potential fixes the magnetic charge directly. 

By integrating the Killing identity V a VbK a = Rb c K c , for the Killing field K a — 8f (which gives R\ = 
— (\ / —gg tt g TT gtt,T)'/\/—g)), together with the Einstein equation R\ = (R — 2A)/2 — 8ttGT}, it is possible 
to isolate the bulk action contribution at r — > oo and on the cosmological event horizon. Here we use the 
observation that, for monopole solutions (A t — 0), the term T\ exactly cancels the matter field lagrangean 
in the bulk action Cm- The integration is from the cosmological horizon out to some fixed r that will be 
sent to infinity. Therefore we shall work in the "upper patch" outside of the cosmological horizon. The 
divergent contribution given by the surface integral term at infinity in R\ is also cancelled by /surface + let 
and we arrive at the simple finite expression 

U = f (M + r 3 c /f). (63) 

2 However, we expect all AF configurations to present dS counterparts, in particular the axially symmetric monopole- 
antimonopole solutions, with a vanishing net magnetic charge Q m = 1411 . Such solutions exist also for A < I3UI . 
3 This definition has nothing to do with the finiteness of the mass-energy of monopole solutions. 
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Since d/dt is a Killing vector, we may periodically identify it and assign a value for p. If we analytically 
continue t — ► it, we obtain a metric of signature (—2,2). The submanifold of signature (— ,— ) described 
by the (t, r) coordinates will have a conical singularity at r = r + unless the i-coordinate is periodically 
identified with period 



0H = 



47T 






-1 


{-a 2 (r)N'{r)) 


r—r c 


4tt/2(t) 


r—r 



(64) 



This is the analogue of the Hawking temperature outside of the cosmological horizon. The relation S = A c /A 
results straightforwardly from (|54|l . 

The overall picture for black hole solutions is, however, much more complicated than in the particle-like 
case. This is because both black hole and dS space produce thermal radiation due to quantum effects. If 
these temperatures are different (which is the generic case, see Figure 6), then the energy flows from the 
hotter horizon to the cooler one and the black hole will gain or lose mass. 

In deriving these relations we supposed that the mass function m(r) approaches a finite value as r — > oo. 
This is the case for the monopole solutions in the EYM theory. However, we could not find EYMH solutions 
with this feature and they are unlikely to exist. For EYMH configurations, the constant M in l|63|) is formally 
infinite, which leads to a divergent action, too. 

For AdS solutions, in some cases it is still possible to obtain a finite mass by allowing the regularising 
counterterms to depend not only on the boundary metric but also on the matter fields on the boundary 
[43U44) . We expect that this approach would apply also to dS EYMH monopoles, yielding a finite mass/action 
to these solutions (see e.g. the case of dS Goldstone solutions 

Also, as remarked in |46| . while the counterterms are necessary to render both the action and energy 
finite, the thermodynamics is not affected by changing the counterterm prescription (basically both I c i and 
(3VJI are shifted with the same amount). Therefore we expect the entropy of the EYMH monopoles, evaluated 
outside the cosmological horizon to satisfy also the generic relation S = A c /4. 

3.3 Remarks on the boundary CFT 

As in the AdS/CFT correspondence, the metric of the manifold on which the putative dual Euclidean CFT 
resides is defined by 

i 2 

7y = lim —//,,. 

The dual field theory's stress tensor t£, is related to the boundary stress tensor by the rescaling |¥7] 

Vrf^ = lim Vhh ik T jk . 

T — >00 

For the bulk metric ansatz (jJJ) , the geometry of the manifold on which the dual CFT resides is given by the 
cylinder metric 



ds 2 = - flj dx l dx j = dt 2 + £ 2 dn 2 . (65) 



The CFT's stress tensor is 



1 



T. 



, M(-3u i « i +$. (66) 



8ttG 

where u l = 8\. This tensor is finite, covariantly conserved and manifestly traceless, as expected from the 
dS/CFT correspondence, since even dimensional bulk theories are dual to odd dimensional CFTs which have 
a vanishing trace anomaly. 

One can see that a divegent 971 in the bulk would imply a divergent value for the energy of the dual CFT. 
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3.4 Quantities inside cosmological horizon 



Following UJ, one may also define a total mass M c inside the cosmological horizon. This can be done by 
integrating the Killing identity V^V^ATp = R vp K p , for the Killing field K = d/dt on a spacelike hypersurface 
£ from the origin to r c to get the Smarr-type formula 

M C =^J V^dE"" = ^qJ v AJf M dE" + f (2T^ - Tg^KOdE?, (67) 

where one integrates between the origin and the cosmological horizon. It is natural to identify the left-hand 
side as the total mass within the cosmological horizon. The first term on the r.h.s. is the contribution of 
A to the total mass within the cosmological constant, while the last one is the contribution of the matter 
fields. For the metric ansatz Q, M c can also be rewritten as M c = — k c A c /AttG = — r^<r(r c )N' (r c ) /2G, 
where k c , A c are the cosmological horizon surface gravity and area, respectively. 
A relation similar to i|67[l can also be written for black hole solutions 

M C S ^ / V„tf ^ = ^ + ^cj Vi + (2T„ - T^)K^. (68) 

One should notice that M c in the above relations is different from 9JI, since the former keeps track only of 
the mass within the cosmological horizon, whereas the latter has contributions also from outside. 

Also, the EYMH monopole solutions we found have always a finite M c , while DJl generically diverges. 

4 Numerical results 

Although an analytic or approximate solution appears to be intractable, we present in this section numerical 
arguments that the known EYMH solutions can be extended to include a positive cosmological constant. 
Apart from that, we discuss also monopole solutions of the EYM model, emphasizing their properties at 
timelike infinity. 

To integrate the field equations, we used in all cases the differential equation solver COLSYS which 
involves a Newton-Raphson method 48 . The case A > leads to the occurrence of a cosmological horizon 
at r — r c with N(r c ) = 0. Similar to the U(l) case, the cosmological horizon radius r c is a function of A. 
Numerically, we have implemented a supplementary equation for the "function" A with A' = 0. Solving 
the set of equations fixing r c then gives us the solutions for a numerically determined value A = const.. 
Using this trick, we can impose two further boundary conditions for the "function" A, one being effectively a 
condition on m at r c . Our procedure then has been to first solve the equations on the interval [0 : r c ] fixing 
r c and then solve the equations on the interval [r c : oo). In all our numerical computations it then turned 
out that the solutions on [0 : r c ] and [r c : oo) can be combined to a continuous solution on [0 : oo). 

4.1 Magnetic monopole solutions in EYMH theory 

The system of equations depends in this case on three parameters (A, a, (3) in the case of gravitating 
monopoles and additionally on r/, in the case of non-abelian black holes. 

Using the initial conditions at the origin/event horizon and on the cosmological horizon, the equations of 
motion (f X3|> - (|14f> were solved varying A for a range of values of the coupling parameter a and several values 
of /3. 

We here present our results for [3 — 0.1 and remark that the results are qualitatively similar for other 
values of (3 we considered. 

As a general feature of both particle-like and black hole solutions, we note that while a negative cos- 
mological constant exerts an additional pressure on configurations, causing their typical radius to become 
thinner 0|, a positive A has the opposite effect, causing the typical solution radius to expand beyond the 
value it would have in asymptotically flat space. 

Hence, for small values of A, the solutions in the region r < r c resemble the AF solutions, being surrounded 
by a cosmological horizon and approach dS geometry in the asymptotic region. 
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Also, as a increases, the cosmological horizon shrinks in size. The dS configurations are generally not 
confined inside the cosmological horizon, with all variables and their first derivatives extending smoothly 
through the cosmological horizon. 

Note that the bounds for finite energy solutions Mjj < Ms and Myy < Mb after the rescalings read 
2/3 < \/3A and y/12 < y/I. 

4.1.1 Solutions with a regular origin 

When A is increased from zero, while keeping a, fixed, a branch of dS solutions emerges from the corre- 
sponding asymptotically flat configurations. This branch ends at a maximal value A max . The value of A max 
depends on the parameters a, (3. For example, for solutions with /3 = 0.1 in a fixed dS background (a = 0), 
we find A max 0.035. The value of A max varies only little with a and /3, e.g. for a = 1, (3 = 0.1 we find 
Amax ~ 0.033. The profile of a typical first-branch solution is presented in Figure^a). 

A second branch of solutions always appears at A max (such that at A max the two branches merge at 
one solution), extending backwards in A to a zero value of the cosmological constant. The trivial solution 
(f)(r < oo) = 0, w(r < oo) = 1 with (f>(r = oo) = 1, w(r = oo) = is approached (the solution thus has an 
infinite slope at r — oo) in this limit. This can be understood as follows: since we have rescaled A — > A/e 2 r/ 2 , 
the limit A — * on the first branch corresponds to asymptotically flat space with r\ ^ 0. On the second 
branch, however, the limit A — > corresponds to asymptotically de Sitter space with r\ — > 0, which then 
implies the trivial solution. We give the profile of a typical second-branch solution in Figure ^b). 

The main difference between the solutions on the two branches is that the matter field functions attain 
their asymptotic values much quicker for the solutions on the first branch. 

This is shown in Figure El where some numerical data is given as function of A for the two branches. 
The results in this figure are obtained for a = 1, j3 = 0.1 but they remain qualitatively the same for all 
gravitating solutions we have considered. The maximal values of A are always below the critical value found 
for solutions in fixed dS background, and as a result the mass of our solutions measured at future/past 
infinity always diverges. While we have found solutions with 2/3 < V3A (finite contribution from the Higgs 
field) , all solutions have v'A < \/l2- Thus although the contribution from the Higgs field might be finite, the 
contribution from the gauge field, however, leads always to an infinite mass. For small values of A (Figure 
la), the divergence in m(r) will manifest itself only for very large values of r; however, this becomes obvious 
for large enough values of A (see Figure lb). 

Note, though, that the mass M c within the cosmological horizon stays finite (see also Figure and that 
the mass M c doesn't tend to zero in the limit of trivial solutions. 

The existence of other disconnected branches of solutions for A > A max appears unlikely. 

For a fixed A < A max , our numerical analysis shows that regular solutions exist on a finite interval of 
the parameter a (which depends on /3), i.e. a <E [0, a max (A)}. Correspondingly, the cosmological horizon 
depends slightly on a : r c (a = a max ) < r c < r c (a = 0). For e.g. A — 0.003, we find 43.6 < r c < 44.7, as 
shown in Figure [3J 

When the parameter a is increased, the local minimum of the metric function N(r) becomes deeper and 
deeper while at the same time the cosmological horizon decreases slightly. In the limit the local 

minimum of N approaches zero at r = rh with N(r = ru) = and A^'(r)| r=rh = 0. In fact this value is 
determined cxplicitely according to 



These results give strong evidence that in the limit a — > a max (A) the gravitating monopole solutions 
bifurcate with the branch of extremal RNdS solution on the interval [r^,oo]. The gravitating monopoles 
separate in this limit into an interior region r S [0 : rh] with a smooth origin and a nontrivial YM field, and 
an exterior region r £ [rj : oo] where the solution corresponds to an extremal RNdS solution with w = 0, 
(f) = 1 and thus unit charge. As a consequence for r close to and larger than the cosmological horizon r c 
the metric functions cr(r), N(r) are identical to ones of the RNdS solution. The mass value can then 




for —r > A > 

4a 2 
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be determined explicitely by imposing the condition N(r c ) = 0. We also observe that in the limit a — > 
the value cr(0) decreases considerably, while the value r c tends to a finite (but A — depending) value. The 
solution corresponds to an hyperbolic monopole. 

As stressed above, the value of the cosmological horizon depends only slightly on a. This is shown in 
Figure |21 We found also that M c and a(0) depend rather weakly on A. 

4.1.2 Black hole solutions 

Similarly to the globally regular solutions, non-abelian black holes solutions can be constructed with A > 0. 
These solutions correspond to non-abelian black holes sitting inside the core of de Sitter monopoles. They 
have a double horizon: an event horizon at r — r n and a cosmological horizon at r — r c . The function N(r) 
is zero at these two points. A generic de Sitter black hole solution is presented in Figure 0] for a = 0.5, 
A = 0.0006 and = 0.1. The event horizon is fixed by hand while the cosmological horizon depends on the 
choice of a, A. r c is then determined numerically. 

For fixed r^, and sufficiently large a the metric function N(r) develops a local minimum at r — r m i n , 
fh < Tmin < r c and the value N(r m i n ) decreases and tends to zero for a maximal value of a. This is seen 
in Figure El for A = 0.003 and r n = 0.3. Very similarly to the regular case, in this critical limit, we end 
up with a solution presenting three horizons. The intermediate one corresponds a double zero of the metric 
function N(r) with N(r = r m i n ) = 0, N'(r)\ r=rmin = 0. On the interval D" m j n ,oo) the solution approaches 
an extremal RNdS solution, while on the interval [r n , r min ] the non-abelian black hole solution remains non 
trivial. 

This phenomenon is illustrated in Figure [5] for A = 0.003 and r n = 0.3. 

At each horizon of these black holes solutions, the surface gravity (determining the entropy) can be 
computed. In terms of the spherically symmetric ansatz, it is given by 



The various quantities characterising the black hole are shown in Figure |SJ In particular, we see that the 
surface gravity at the inner event horizon decreases considerably when increasing a and becomes very small 
(along with a(r n )) in the limit a — > a max . In contrast, the surface gravity at the cosmic horizon varies only 
little with a (along with u{r c ) = 1). 

Similar to the case of solutions with regular origin, the mass function m(r) of the black hole solutions 
also diverges at infinity, while the mass inside the cosmological horizon stays finite. In order to demonstrate 
the peculiar behaviour of the Higgs function 4>(r) for r — > oo, we have to construct solutions for large values 
of A. In Figure 7, we show our results for A = 0.033. We especially show the value of rip' /(l — <fi). The 
asymptotic value of this particular combination (which tends to some constant ~ 0.75 for the parameters 
chosen here) confirms clearly the behaviour (|32J) . Of course, this phenomenon is present for generic values 
of the parameters although not always easy to set up numerically. We were able to exhibit this property for 
regular solutions as well. 

4.2 Magnetic monopoles in EYM theory 

We have integrated the system of eqs. (|15H for a range of A by using the same techniques. Both solutions 
with regular origin and cosmological black holes have been considered. 

All numerical solutions we found have w 2 < 1, which implies a nonzero node number of the gauge 
potential w. This can be proven by using the sum rule |21| 



which follows directly from the YM equations (ro here is tq — for solutions with regular origin or = r n , 
for black holes). Suppose that w(r) never vanishes and w 2 < 1 for r < r < r c . Then l.h.s. of the above 
relation vanishes, while the integrand of the r.h.s. is positive definite. Therefore the gauge potential of the 





(69) 
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nontrivial YM configurations with w 2 < 1 must vanish at least once in the region inside the cosmological 
horizon. 

Thus for a given value of A, the solutions are indexed by the node number k. Different from the AF case, 
however, the asymptotic value of the gauge potential is not fixed by finite energy considerations. Similar to 
M, Wq appears as a result of the numerical integration. 

For small values of the cosmological constant A -C 1/-R 2 {R c corresponding to the typical solutions' 
radius), the contribution Ar 2 to the energy density is negligible. For r <C R c , the solutions do not considerably 
deviate from the corresponding flat space configurations. In the region r > R Cl the effects of A become 
significant. 

Several characteristic features of the two node particle-like solutions of the EYM-A system are plotted as 
a function of the cosmological constant in Figure 8. One can see that both the mass inside the cosmological 
horizon and the mass evaluted at timelike infinity are finite. The next two plots show the profiles of typical 
EYM-A solutions with regular origin (Figure 9) and cosmological black holes (Figure 10). In both cases, 
u> 2 (oo) ^ 1, which implies the existence of a nonvanishing nonabelian magnetic charge defined according to 

One should remarks that the EYM theory presents solutions with dS asymptotics only for values of the 
cosmological constant up to some A max < 3/4 22 . 

As discussed in , the dS EYM-SU(2) theory is a consistent truncation of D = 11 supergravity via 
Kaluza-Klein dimensional reduction on a non-compact space, the positive cosmological constant being fixed 
by the SU(2) gauge coupling constant A = 4e 2 . The "internal" space is a smooth hyperbolic seven-space 
written as a foliation of two three-spheres, on which the SU(2) Yang-Mills fields reside. However, no solutions 
with dS asymptotics exist for this value of the cosmological constant 



5 The nonabelian monopole in cosmological coordinates 

Although the line element J2J takes a simple form in a static coordinate system, the expansion of the universe 
is not manifest. Furthermore, the static coordinates (t,r) break down at the Killing horizons. However, one 
can prove that the results in the previous sections remain qualitatively unchanged when using a different 
parameterisation of dS spacetime. A particularly interesting case is provided by the planar coordinates (for 
which the spatial slices are flat) — or cosmological or inflationary coordinates — with a dS line-element 

dg 2 _ e 2HT^ dR 2 + R 2 dQ 2^ _ ^ 

(to conform with the standard notation in literature, we note here A = 3H 2 , i.e. H = l/£). The relation 
between the physics in a static coordinate system and in a planar coordinate system in the case of Einstein- 
U(l) theory is discussed in [3~3"| . 

Considering the general metric ansatz (|7|) and the matter fields form @, (|10f> . the solution in the 
cosmological form can be written as: 

ds 2 = -V- 2 dT 2 + U 2 e 2HT {dR 2 + R 2 dfl 2 ), 

' w(i?e HT )rid6»+ (cos 9t 3 + uj(Re m )T 2 sin 6)dip 



A 

2e 



$ = (j)(Re Ht )T 3 (70) 



where U and V are functions of p = Re HT . The coordinate transformation that relates to the above 
metric form is: 

U( P ) = T M, t = T -[ Hrdr (71) 

P J a{r)N{r)yf<j{r) 2 N{r)+H 2 r 2 

V- 2 (p) = a 2 {r{p))(^±^-\\ (72) 



while r(p) is determined implicitly by: 



U 



a(r)dr 



r^fa 2 {r)N{r) +r 2 H 2 



= lnp. (73) 
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As an example, let us consider the U(l) solution with: 



w(r) = 0, a(r) = 1, 0(r) = 1, 
Then we can easily integrate (I7r>|) and we obtain: 

M 2 -a 2 



N(r) = 1 + — - H 



U{p) 



1 + H 
P 



Ap 2 



1 + M 



M 2 



4/< 2 



1 



Af 2 — Q 
4p 2 



(74) 



(75) 



which is readily seen to be the RNdS solution in cosmological coordinates. The extremality condition can 
be written in this case as |a| = M. 

One advantage of using the planar coordinates is that the expansion of the universe is now manifest: 
the cosmological expansion chart corresponds to H > 0, while the contracting chart is given by H < 0. 
However, this comes with the price that the manifest time-translation symmetry is broken while the charts 
that cover the horizons are highly distorted. The boundary geometry Z ± is now approached for large T and 
the boundary topology is R 3 (written above in spherical coordinates). 

Let us notice that p — Re HT is unchanged under the transformations: 



T 



T 



R 



a R, 



(76) 



which means that the geometry in the cosmological ansatz is preserved. The first term in the above transfor- 
mation generates the time translations in the bulk, while the second term corresponds to a scale transforma- 
tion of the radial coordinate on the boundary. There exists a Killing vector associated with this symmetry 
and it can be written as: 



(77) 



Similar to the (electro-)vacuum case, the norm of this Killing vector will vanish precisely where N(r) = 0, 
that is, at the horizons r — r/j, r c as determined in static coordinates. 

It is straightforward to compute the boundary stress-tensor in cosmological coordinates. Choosing the 
boundary dAi to be the three-surface given by a fixed value of T, the normal to this surface is = I/VSt^ 
we find the following components of the extrinsic curvature: 



K 



HU 2 a 2 



RR — 



K, 



HU 2 a 2 r 2 



oo 



K, 



HU 2 a 2 r 2 sin 2 9 



iptp 



(78) 



that is, Kij = ——hij, where hij is the induced metric on dM. and we denoted a(T) 



The components of the boundary stress-tensor in cosmological coordinates are then given by: 



SttGTrr — 



87rGTfl ; 



2U 2 a 2 (H 
2U 2 a 2 R 2 



1 

~ I 
H 1 



8ttGT^ 



a" \ a 
2U 2 a 2 R 2 sin 2 6> 



H 
a 



2£U h 
rU 

-ir 1 

1 

7 



RU R \ 


2Mi 


2U J 


R 3 a 


¥f- 


U.R 




h RU + 



£R 2 sin 2 



U 



+ 0(a- 2 ) 

U.rr 
U 

Ur 
+ — 

rU 



2M£ 
Ra ' 

Urr 



U 



0(a- 2 ) 

2M£sm 2 9 



Ra 



+ 0(a~ 2 ) 



where we assumed an expansion of the form U(p) = 1 + + O(^). Notice that in the limit T — > oo the 
first terms on the right-hand side should cancel out once we take advantage of the fact that in this limit 
a -> 1 (recall that H =1/1). 

The conserved quantity associated with the Killing vector £ will be the mass and it can be evaluated at 
past or future infinity according to the sign of H. In our case the conserved mass is readily seen to be 



Tl = 



-M, 



(79) 
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which coincides, as expected with the value in a static coordinate system. In particular, a divergent 9Jt is 
found again for Mjj > Ms or Mw > Mb- 

According to the dS/CFT conjecture, the dual field theory should live in our case on an Euclidian 
manifold whose metric is : 

ds 2 = li] dx l dx j =<1R 2 + R 2 dn 2 (80) 

obtained by an infinite conformal rescaling of the boundary geometry on dA4 , corresponding to a flat D = 3 
Euclidean metric. The stress energy tensor of the dual theory will be given by: 

^ G -R = ^ 8^ = 8^ = -^. (81) 
It is easy to see that this stress-tensor is covariantly conserved and traceless, as expected on general grounds. 



6 Conclusions 

This work was partially motivated by the question of how a positive cosmological constant will affect the 
properties of a gravitating monopole. To the best of our knowledge, this question has not yet been ad- 
dressed in the literature, except for the preliminary results in Apart from this motivation, the study of 
gravitating matter fields configurations in asymptotically dS space may help a better understanding of the 
conjectured dS/CFT correspondence. 

We have found that despite the existence of a number of similarities to the A = case the asymptotically 
dS solutions exhibits some new qualitative features. All solutions present a cosmological horizon, on which 
matter fields take non-trivial values. An interesting feature of the dS solutions appears to be the absence of 
monopole configurations without a Higgs potential. Also, contrary to the naive expectation that a small A 
will not affect the properties of the configurations drastically, we find that the mass of dS solutions evaluated 
at future/past infinity by using the quasilocal tensor of Brown and York diverges (although the mass within 
the cosmological horizon stays finite). 

Since this result is based only on the asymptotic expansion of the field equations together with the Higgs 
mechanism, we expect it to remain valid for the dS versions of various possible extensions of the magnetic 
monopole model. Moreover, our results appear to be a generic feature of the asymptotically dS particle-like 
solutions of a spontaneously broken gauge theory. Similar qualitative results are found for dS sphalerons 
with a Higgs doublet |T3). 

The question of the stability of our solutions is very important. For the monopoles, we know that the 
solutions on the first branch are stable for A = 0, since these are nothing else but gravitating monopoles 
in asymptotically flat space. Invoking arguments from catastrophe theory [20], we thus believe that the 
monopoles on the first branch for < A < A max are stable, while the solutions on the second branch have 
one mode of instability. This statement should definitely by checked by means of a detailed analysis, e.g. 
with the one discussed in [ST] , 

Supplementing the EYM system with a cosmological constant leads to the occurrence of the "bag of 
gold" family of solutions, where an equator naturally appears in the metric and the spacetime becomes 
spatially compact and homeomorphic to a three-sphere [22] ■ Nothing like that was observed in the present 
context because the Higgs function is imposed to approach its expectation value asymptotically. Relaxing 
this condition, however, we think that it could lead to compact solutions as well for fine tuned values of the 
cosmological constant like e.g. in We plan to study this possibility for EYMH in near future. 

Also, it would be interesting to apply the isolated horizon formalism, recently extended to asymptotically 
dS spacetimes to the particular system discussed here. 

The divergencies described above are not too disturbing for physics inside the cosmological horizon. 
However, a divergent mass-energy and action appear to lead to severe problems concerning the possible 
holographic description [H] of a gravitating spontaneously broken gauge theory in dS spacetime. A divergent 
ADM mass has been found also for solutions of some theories in asymptotically AdS spacetime. However, 
in some cases it is still possible to obtain a finite mass by allowing the regularising counterterms to depend 
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not only on the boundary metric but also on the matter fields on the boundary |44| . It would be interesting 
to generalise this method to the dS case and to assign a finite mass (evaluated outside the cosmological 
horizon) to the solutions of a spontaneously broken gauge theory. We believe that this may lead to further 
understanding of the rich structure of a field theory in dS space as well as profound implications to the 
evolution of the early universe. 

It is likely that a re-examination of various field theory nonperturbative effects for a dS ground state may 
lead to further surprises. 
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Appendix A: BPS solutions with a Liouville potential 

In it was shown that the BPS monopoles of four-dimensional YMH theory continue to be exact so- 
lutions of this model even after the inclusion of gravitational, electromagnetic and dilatonic interactions, 
provided that certain non-minimal interactions are included. The solutions presented in |27| are asymptot- 
ically flat. Here we prove that a similar construction can be done for a Liouville dilaton potential, i.e. a 
positive effective cosmological constant. 

Following the conventions in |27) . we consider a four dimensional action principle on the form 



S 



4ttG 



R 1 



(W) 



P 2ha Lf tv v 



-Ibo 



-A 



Smat 



(82) 



with a the dilaton field 4 b the dilaton coupling constant and / M „ the U(l) field. Different from the situation 
in [22], the above action principle contains a Liouville potential V(a) = — e~ 2bCT A/4. The expression of the 
matter action is given by 



Smaller = jd X\/ Zr g\ - -C 



1 (l-i> 2 ) I <l+i> 2 ) Cvl+6 2 a I 

" a Tr(F^F» v ) - -e-^ a Tr(D^D^) + ^= f^F a0 ) \, (83) 



with F^y the SU(2) field, $ the Higgs scalar, while c is a constant defined bellow. Also, to simplify the 
general picture, we take here the gauge coupling constant g = 1. 
The Einstein-Maxwell-dilaton field equations are given by 



G>„ - 2T^(/) - 2T M „(<7) - (8wG)T^(mat.) 

V JU (e 2bCT /^) = (4nrCf)J v (mat.) 



v y L J J [iv 



(84) 



where we note 



V(mat) = e ^ a Tr ( F^\F V - -g^F a/3 F 



a/3 



T Tr ( D^Dv® ~ -g^D a <S>D a <^ 



J v (mat.) 



Vl + bV^tr^G^) 



(85) 



'The dilaton a should not be confused with the metric function cr(r) which enters the ansatz 171 , 
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In 127] it has been shown that for A = the equations of motion for and <& are solved by any YMH 
configuration that solves the fiat space Bogomol'nyi equations, 

Gij = T e^D^ , (86) 

provided c = 1. To this aim, a special metric ansatz has been used, in terms of only one metric function. 
The the Einstein, Maxwell and dilaton equations are then equivalent to a single Euclidean 3-space Poisson 
equation. In the absence of matter fields S m atter — 0, this choice gives the "extreme" dilaton black hole 
solutions of 

This construction can be generalised for a Liouville potential term in the action principle. However, given 
the presence of an effective cosmological constant, the metric ansatz would be time-dependent. 
A straightforward generalisation of the metric ansatz used in |27| is 

ds 2 = -e- 24, dt 2 +R 2 (t)e 24 'dx-dx, where = e a ' b /R{t), R(t) = {-) 1/b \ (87) 

to 

the only nonvanishing component of the U(l) potential being 



VT+¥ y ' 



The free parameter to is fixed by A 



It is convenient to define 

e^ = (l + i?r( 1 + h2 )/(x)) 1 /( 1+b2 ). (90) 

A direct computation shows that, for c 2 = 1, the YMH equations implied by the action principle (|83[1 are 
automatically satisfied by any flat space solution of the Bogomol'nyi equations. 
Then, the function /(x) is uniquely determined as a solution of the equation 

V 2 / = -47rG(l + 6 2 )tr(A$-A$), (91) 

which follows from the Einstein-Maxwell-dilaton equations. In the absence of S mat teri one recovers the 
(multi-)black hole solutions in cosmological Einstein-Maxwell-dilaton theory discussed in |54| . 
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Figure 1: The profiles of two typical gravitating EYMH monopole solutions with a regular origin are plotted 
as a function of r for a = 1. In (a), we show a solution on the first branch for a = 0.5, in (b) a solution on 
the second branch for a = 1.0. Here and in figures 171 1§1 and UHl we have indicated also the position of the 
cosmological event horizon. 25 



1.5 
1 


I 

a=l (3=0.1 




i 

1 


1 


i i 

(a) 


0.5 


- 4>'(0) 


i 













0.5 








1 




-1 








... -f(r c )xl0 2 


1.5 

O 






2 

1 







0.005 0.01 0.015 0.02 0.025 0.03 0.035 

A 





A 



Figure 2: The dependence of solutions' properties on the value of the cosmological constant is plotted for 
particle-like monopole solutions. The labels 1 (rcsp. 2) refer to the first (resp. second) branch. 
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Figure 3: The mass inside the cosmological horizon M c , the value of the metric function a at the origin, 
cr(0), and the value of the cosmological horizon r c of the gravitating monopoles are given as functions of a 
for fundamental branch solutions with A = 0.003. 
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Figure 4: The profiles of the functions of the non-abelian black holes are given for a = 0.5, r/j = 0.1 
A = 0.0006. 
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Figure 5: The metric functions TV and a of two EYMH black hole solutions with rt — 0.3, A = 0.003 and 
two different values of a are shown as a function of r. 
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Figure 6: The surface gravity at the two horizons, the mass inside the cosmological horizon M c and the 
minimal value N m of the metric function TV, are shown as functions of a for non-abelian EYMH black hole 
solutions with r h = 0.3, A = 0.003. 
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Figure 7: A typical black hole solution in EYMH theory is plotted for a = 1, ru = 0.3 and A = 0.033. We 
also plot rcf>' /(l — 4>') which indicates the power with which <f> reaches its asymptotic value. 
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Figure 8: Several characteristic features of the two node particle-like solutions of the EYM-A system are 
plotted as a function of the cosmological constant. Here M and M c are the values of the mass function m(r) 
at infinity and at the cosmological horizon, r c is the event horizon radius and tr(0) is the value of the metric 
function a at the origin. Also, wo and w(r c ) are the values of the magnetic gauge potential at infinity and 
at the cosmological horizon. 
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Figure 9: A typical asymptotically de Sitter monopole solution with a regular origin in EYM-A theory. The 
value at infinity of the magnetic gauge potential is wq — 0.761. 
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Figure 10: A typical asymptotically de Sitter black hole monopole solution in EYM-A theory The cosmo- 
logical horizon radius is r c = 15. The value of the gauge field on the event horizon is Wh = 0.9936, while at 
infinity w = 0.8244. 
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